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1. Introduction 



Rare B decays proceeding through b — > s flavour changing neutral current transitions provide 
particularly sensitive probes for physics beyond the Standard Model. Experimental uncertainties 
in B — > K*y have been reduced to the few -percent level, and those for B — > K^Hi are presently 
around 10 percent [p. So far, measurements are consistent with Standard Model expectations 
[0, [|]. Future experiments like LHCb are going to reduce the experimental uncertainties further. 
Thus, we are entering an era of precision flavour physics, requiring experimentalists and theorists 
to focus on reducing statistical errors and employing multiple cross-checks to further quantify 
systematic uncertainties. 

Lattice QCD can contribute to this by providing first-principles nonperturbative calculations 
of hadronic matrix elements relevant for rare B decays. In this report, we present our current 
progress in doing these computations. We use a non-relativistic effective action for the b quark 
called moving NRQCD (mNRQCD) [Q], which helps reduce discretisation errors in the final state 
meson at large recoil. This is achieved by giving the B meson a significant velocity in the lattice 
frame, so that at a given q 2 the momentum of the final meson is smaller. 

Previous progress of this project was reported in In this work, we investigate how far 
statistical errors can be reduced with the help of stochastic (random wall) sources, and exam- 
ine different fitting strategies. Furthermore, all calculations are now done with the full &{\/m 2 b ) 
mNRQCD action exactly as in [Q]. 

This report is organised as follows. We discuss the phenomenology of rare B decays and 
our strategy for the lattice calculation in Sec. |2| Details of the operator matching for heavy-light 
currents with mNRQCD are given in Sec. ||. Then, in Sec. ^ we explain the use of stochastic 
sources and compare their effectiveness to point sources using numerical calculations. Finally we 
show some preliminary results for the form factors in Sec || and conclude in Sec. 



2. Physics of rare B decays and strategy for the lattice calculation 

The starting point for studying weak decays of hadrons is the effective weak Hamiltonian. For 
b — > s transitions in the Standard Model, the governing Hamiltonian is 1 

= -^v^EqoogKm) (2-D 

where Gp is the Fermi constant, C;(jU) are Wilson coefficients taking into account short distance 
physics and Qi are effective local operators. In the Standard Model, there are 10 operators we need 
to consider for radiative and semileptonic decays: 

Qi = (siCj) V - A (cjbi)v-A Qi = {sc) V - A {cb) V -A 

Q 3 = {sb) V -A"L q {qq)v-A Qa = {sibj) V -A"L q {qjqi)v-A 

Q 5 = (sb) V -ALq(qq)v+A Qe = {sibj)v-Al< q {qjqi)v+A 

Qi = ^mbSiG^il+^biF^ Q 8 = JL imb s i c^{l + %)T i a j bjGl v 
Q 9 = &(sb)v- A (te)v Gio = &(sb)v-A(Ze)A 



To obtain Eq.(|2.1|), V u bV* ls <C V^V,* and CKM unitarity are used. 
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Figure 1: Penguin and box diagrams governing b — > sy and b — > stt in the Standard Model. (Decay of 
photon or Z to It not shown in the penguin diagram.) 



Figure 2: Decays b — > 57 and ft — > s££ via short distance operators in the effective weak theory. 




■ B \ 




Figure 3: Long distance contributions can come from charmonium resonances via Q2 (left) or weak annihi- 
lation (right). (Weak annihilation is doubly CKM suppressed fori? — ► K*.) 



The dominant short-distance contributions come from Q-j, Q<j and Q 10. They arise from penguin 
and box diagrams illustrated in Fig. |]. The hadronic matrix elements of these quark bilinears 
(Fig. ||) are straightforward to compute in lattice QCD, at least for some values of initial and 
final state momenta. The matrix elements of the quark bilinear currents in Qj, Q9 and Q\q are 
parametrised by form factors as follows: 



(P(p')\qfb\B(p)) = f + (q 2 



, f ( „2, M 2 B -Mj 

+Mq ) — - — 9 ) 



qv (P(p')\q0^b\B(p)) 
(V(p',e)\qfb\B(p)) 



2iV(q 2 ) 

Mg+My 



F Hvpa * 1 



(V{p',e)\qf%b\B{p)) = 2M v A (q 2 )^-q^ + (M B +M v )A l (q 2 ) 



2 e-q 



p^+p 



'11 



M 2 B -M 2 
q 2 q 



M B + M V 

q v (V(p',e)\qa^ v b\B(p)) = 4T l (q 2 )e^ pKa e* p p K p' (y , 
q v (V{p',e)\q^ v f 5 b\B{p)) = 2iT 2 (q 2 ) [e* (M 2 -M 2 ) - (e* • q)(p + p\ 



+2iT 3 (q 2 )(e*-q) 



111 



M 2 B -Ml 



(p + p'h 



(2.2) 
(2.3) 
(2.4) 

(2.5) 
(2.6) 

(2.7) 
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where e is the polarisation of the vector meson and q = p — p'. 

The form factors are functions of q 2 . We would like to compute them for the whole range 
of q 2 directly using LQCD; however it is not possible to do so on current lattice ensembles. Our 
calculations are most reliable in the low recoil limit q 2 <?max> where both the B meson and the 
final state meson are roughly at rest. The march to high recoil (small q 2 ) is blocked by three 
barriers: discretisation errors as spatial momenta in the lattice frame become comparable to the 
inverse lattice spacing, errors in heavy quark effective theory as the QCD dynamics of the decay 
become comparable to mt, rather than Aqcd, and growing statistical errors. 

Our strategy for calculating the q 2 -dependence of the form factors is the following: (1) Com- 
pute with large values of q 2 where the form factors above are the dominant hadronic contributions 
to B — > K(*>£1 decays. The £ = /I decays are clean channels at the LHC, and Standard Model cal- 
culations of the differential cross section with q > m 2 c will be tested directly against experimental 
measurements. (2) Use a combination of methods, including moving NRQCD, to push LQCD 
calculations to smaller values of q 2 , determining the shape of the form factors over a range large 
enough to test various phenomenological ansatze for the shapes. (3) Use these fits to extrapolate to 
the low q 2 region, down to q 2 = 0. 

We have discussed above the contributions from Qj, £2g, and £2io- Compared to these, con- 
tributions from Q\ and £23-6 are loop-suppressed. The contribution from Q& is also suppressed 
by an approximate factor (Aqcd/wibXQ/C?) ~ 0.05. There are, however, complications which 
have to be understood due to £22- The contribution from £22 dominates when the cc pair is resonant 
(Fig. ||, left). This is non-local and we can not calculate it on the lattice. Therefore, we will not be 
able to extract the differential branching fraction in ranges of q 2 which are close to a charmonium 
resonance. Nevertheless, the form factors we do calculate will be unaffected by these resonances 
and extrapolations in q 2 will be smooth. 

Weak annihilation contributions, where the valence quarks of the B meson annihilate into a W 
which decays into a u quark and a d or s quark (Fig. [| right), are long-distance contributions to 
radiative and semileptonic B — > decays, since the photon or Z is emitted separately from the 
flavour-changing interaction. This process is highly suppressed for b — > s decay since V U },V* S -C 
VtbV t *, in contrast to b — > d decay where there is no CKM suppression. 

Another complication facing LQCD calculations with vector mesons in the final state is their 
instability to strong decays. With the presently used value for the light quark mass, the p and K* 
are stable, but eventually we need to extrapolate to the physical light quark mass. The effects at the 
decay threshold can be perhaps be studied in separate work by looking at the vector meson decay 
constant as the quark mass is reduced. We can also try to check the validity of the extrapolation by 
comparing \V u b\ obtained from B — > plv with that from B — > %tv. 

3. Moving NRQCD and matching of heavy-light operators 

To reduce discretisation errors for the light meson in the final state at high recoil, we work in 
a reference frame where the B meson is not at rest, so that for a given value of q 2 the momentum 
of the light meson is reduced. Moving nonrelativistic QCD (mNRQCD) allows us to treat the 
momentum of the heavy quark arising from the frame choice exactly. 
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Figure 4: Matching coefficients for the vector cur- 
rent. 




frame velocity v 



Figure 5: Matching coefficients for the tensor cur- 
rent. The renormalisation scale is pi = m b . 



The full &{AQ CD /tnfy mNRQCD action used in our calculation was derived in Ref. [Q]. 
Schematically, the 4-momentum of the b quark is parametrised as p = m b u + k, where u = 7(1, v) 
and 7=1/ Vl — v 2 . Then a non-relativistic expansion in the 3-momentum k is performed. 

At tree-level, the QCD heavy quark field f(jc) is related to the mNRQCD two-component 
quark- and antiquark fields y v {x), £ v (jc) by the transformation 



W{x)=S{A)T FWT e 



with 



%{x) 




(3.1) 



where e -" nux f removes the additive heavy quark mass term, Tpwr is the Foldy-Wouthuysen-Tani 
transformation, S(A) is the spinorial representation of the boost and Ap, removes time derivatives 
in the Hamiltonian. 



Eq. (3.1) can be used as a starting point for the construction of the lattice operators correspond- 
ing to the heavy-light currents in the operators Qt, Qg and <2io- For a current of the form / = qFb 
with some Dirac matrix F, one obtains, using the equations of motion to eliminate time derivatives, 



±qrs(A)% 

/(+) 



(+) 



+ 



1 



2m by /y 



q r {-ifv + ij+iv/y) D 5(A) % {+) + ff(l/m 2 b ) (3.2) 



where ¥^ has vanishing lower components. On the lattice, the continuum covariant derivative 
D has to be replaced by a discrete version A. For the light/strange quark q, we use an improved 
staggered action; presently ASQTAD. This means that the field q has to be expressed in terms of 
the staggered field % q . 



Eq. (3.2) does not yet include any radiative corrections. Since the second term in (3.2) is 
suppressed by AgcD/^fo relative to the first term, we currently restrict our calculation of radiative 
corrections to the first term. To obtain these, we note that the spinorial boost is given explicitly by 

1 



5(A)=5 + (A) 



;i + 7)l-7v-77° 



(3.3) 



As can be seen, ( |3.3[ ) contains a sum of two different Dirac structures; these will mix under renor- 
malisation. Thus, at order 0{a s ), we also need to consider the combination with the opposite 
sign: 



S- (A) 



1 



(l+7)i + 7v-f7° 
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The lattice current through order ff(cc s ) then reads 



7 (0)lat = (1 + aiC+)7 (0) + asC _ 7 (0) (34) 



with 

jf = J_qr S+ {A)^ + \ / (0) = ±= q F 5_ (A ) w} + 1 

Vf Vr 



The matching coefficients c± in ( |3.4| ) are obtained by requiring that the matrix elements in the 
lattice theory and in the continuum (MS scheme) are equal to order a s . 

Results for the matching coefficients for the vector (F = y^) and tensor (F = o^ v ) currents are 
shown in Figs. |4] and ||| These were obtained using automated tadpole-improved one-loop lattice 
perturbation theory Here, the bare b quark mass was set to ami, = 2.8 and the boost- velocity 
is pointing in 1-direction. In the figures, the different Lorentz indices of F are indicated as 
(temporal), || (parallel to v) and _L (perpendicular to v). 

The tree-level @{AQCD/ m b) correction is 

/?> = ±J_ r (-if*+#+m-A 5+(A) y v ( + ). (3 .5) 
m b y/y 2 

Because of the mixing-down that occurs, one should work with subtracted 1/m currents in the 
calculation of the form factors. We have already computed the non-perturbative matrix elements of 
(3.5), so that we can include them in the form factor results once we have obtained the necessary 
mixing coefficients from perturbation theory. 



4. Two-point and three-point functions with stochastic sources 

In order to extract the form factors, we need to compute two-point functions for the light meson 
and the B meson, as well as three-point functions with the currents 7 lat (see section ||) inserted. The 
basic method for the extraction of the matrix elements from fits to correlators was outlined in [g]. 

Here, we investigate the use of stochastic (random wall) sources, which allow us to obtain 
approximations to the all-to-all correlators and therefore possibly to reduce statistical errors. For 
the B — > niv semileptonic decay in standard NRQCD (i.e with the B meson at rest), random wall 
sources have been tested in Ref. [^]. We now include vector meson final states and moving NRQCD 
in our study, and also work at smaller light-quark masses. 

For the light valance quarks, we use one-spinor-component staggered fermions. We work with 
four-spinor-component naive quarks in constructing the interpolating fields for heavy-light and 
light-light mesons. The naive quark propagator G q {y,x) is related to the one-spinor-component 
staggered quark propagator G Xq {y,x) by 

G cl (y,x)=G Xq (y,x)n(y)Q^x), (4.1) 

where Q (x) = Ul=o(%Y" ■ 

The numerical calculations described in the following sections were performed using 400 
"coarse" MILC gauge configurations [^] that have V = 20 3 x 64 and aT l rj 1.6 GeV. These 
configurations include 2 + 1 flavours of ASQTAD sea quarks with masses am u = amj = 0.007 
and am s = 0.05. The ASQTAD action was also used for the light valence quarks, with masses 
am u = amd = 0.007 and am s = 0.04. The heavy quark mass and stability parameter were set to 
arrib = 2.8 and n = 2. 
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4.1 Light meson two-point functions 

For flavour non-singlet light pseudoscalar mesons with 3-momentum p , the exact all-to-all 
two-point function (with naive quarks) is given by 



C 55 (t=yo-xo,i?) = ^(0 5 (y)<Pl(x))e- r P'<y-^ 



L? 
1 

L? 

4 
L? 



Note that due to the taste doubling, the amplitude obtained from this correlator is too large by a 
factor of 16 (in the continuum limit), which needs to be divided out. 

Compared to a single-point-source correlator, where one does not sum over x, we would expect 
a reduction in the statistical errors by a factor proportional to ^Jnv> n L? . However, it is forbiddingly 
expensive to compute (O) directly. Instead, one can approximate (4.2) using random wall sources, 
as explained in the following. 

We define G Xq (y,x ,p>) = £ G Xq (y, X ) ^{x)e^ 

X 

where % p (x) is a vector in colour space (colour index not shown explicitly), with every colour 
component at every spatial site an independent Z2 x Zi_ random number (note that in [0] a different 
type of noise was used). The index p labels different random samples on a given configuration; 
additionally it is understood that new random numbers are used on every gauge configuration. The 
noise fields satisfy , n z 

(4.3) 

nz P =i 

where nz is the number of random samples (this relation becomes exact in the limit nz — > °°). 
Therefore, an approximation to the all-to-all pseudoscalar two point function can be obtained as 

Ml0WS: 41* . ^ _ _ _ 

C 5 5,rwM') = 71 - £ L G x j(y,xo,0Y ■G p Xi (y,xo,P>)e- ip - y (4.4) 

L nz p=l y 



In practice, a small number nz is sufficient; even nz = 1 works to improve the signal for pseu- 
doscalars compared to point sources. However, note that for each non-zero value of the momentum 
p', new inversions are required. 

For a flavour non-singlet vector meson with interpolating field <Pj = q'fjq, the exact all-to-all 
two point function is given by 



y,x 



L? 

yj 
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Figure 6: Comparison of relative errors of K meson two point functions from random wall sources (RW) 
and a point source. "RWcd" means random wall sources with colour dilution. See the main text for the 
numbers of inversions used for the different source types. 

where the phase factor of {—\) x > +y i comes from the relation &(x)% fjQ (x) = (—l) X} ¥sfj- 

To obtain the random wall correlator in this case, a factor of (— l) Vj is added to the stochastic 
source for the zero-momentum quark propagator; we define 

fybxoJ) = I G«M (4-6) 

X 

This means that additional inversions are required for the different polarisations j = 1,2,3. The 
random wall correlator is then 



1 «z 

c ;7 ,rwM') = £ - £ IX,()^o,7T • e£(y,W)(-i) w «- <if *. (4.7) 



A modification of the simple random wall sources is colour dilution, where a Kronecker delta in 
colour space is introduced in the stochastic source, i.e. %£(z) is changed to S a ^ ao E, p (z) for each 
source colour ciq. Here, ^ p (z) does not have a colour index. Then for each colour «o one needs to 
perform a separate inversion and in Eqs. ( |4.4| , |4.7| ) a sum over source colour has to be added. 

Fig. ^ shows the comparison of the relative errors of K meson two-point functions from random 
wall sources and a point source at zero momentum and at ap = 2n( 1 , 0, 0)/L. "RWcd" refers to the 
random wall source with colour dilution and "RW" the simple random wall source. The statistical 
errors are significantly reduced by using random wall sources. However the improvement decreases 
as the momentum of the meson increases. We put four sources on each configuration for each 
type of sources. Therefore the number of inversions for the RWcd and point source is the same 
(1600 x 3) while three times fewer inversions are used for the standard RW source. We use the 



Bayesian fitting method described in Ref. [ 1C] to fit the kaon two-point functions with a function 
of the form 



N-l r , N-l 

T55\2 



e -E„t + e -E„(L,-t) 



(4.8) 



n=0 n=0 

where oscillating terms are included for the parity partners. To ensure the correct ordering of the 
states, we actually use the logarithms of the energy differences, \n(E n — E„-i), ln(E n — E n ^i) as 
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Figure 7: Left: A constrained five exponential fit to the K meson correlator from the random wall source 
with colour dilution (fit range: t — 2. ..62). Right: Fit values for the ground state amplitude against the 
number of exponentials N in the fits, for the RWcd and point sources. 



the fit parameters. Also, the amplitudes of the excited states (n > 0) are written as A n = B„Aq, 
A n = B n Ao, and the relative amplitudes B n , B n are used as the fit parameters. The numbers of 
exponentials N and N are increased until the results are stable. For the kaon at zero momentum, 
one can set N = 0. 

An example fit to the K meson two-point function at zero momentum with the RWcd source 
is shown in Fig. [7] (left-hand side). On the right-hand side of Fig. ^, we compare the results 
for the ground state amplitude from the RWcd and point sources. We observe a factor of 4.2 
improvement in the ground state amplitude and 1 .4 in the energy. For the kaon with momentum 
ap = 2^(1,0,0) /L, these improvement factors are 3.6 and 2.6, respectively. There is a small (about 
1 percent or 1.8a) deviation in the amplitudes between the RWcd and point source correlators. Note 
that we inadvertently used different temporal boundary conditions for the two types of sources, but 
we would expect this to have an effect much smaller than 1 percent. 

For the vector meson correlator with zero momentum, we average all three polarisations. With 
nonzero momentum in the i-direction, only the two transverse polarisations are averaged. In Fig. || 
we compare the relative errors of K* two-point functions from random wall sources and a point 
source at zero momentum and at ap = 2^(1,0,0) jh. The number of inversions for the RWcd, RW 
and point sources are 1600 x 3 x 4, 1600 x 4 and 1600 x 3, respectively, for zero momentum. For 
a nonzero momentum, the numbers are 1600 x 3 x 3, 1600 x 3 and 1600 x 3, respectively. 

As can be seen in Fig. |j| for the vector meson two-point functions a reduction in the statistical 
errors is only seen at small t. An example fit for the K* at rest with the RWcd source is shown in 
Fig. H (left-hand side). On the right-hand side of Fig. ^, we compare the results for the ground state 
amplitude from the RWcd and point sources. Unfortunately, no improvement is observed, even 
though the computational cost for the RWcd source was higher. 

4.2 Heavy-light two-point functions 

For B/B s mesons, the exact all-to-all correlator with moving NRQCD and staggered actions is 
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Figure 8: Comparison of relative errors of K* meson two point functions from random wall and point 
sources. See the main text for the numbers of inversions used for the different source types. 




t Number of exponentials N 



Figure 9: Left: A constrained five exponential fit to the K* meson correlator from the random wall source 
with colour dilution (fit range: t — 1...63). Right: Fit values for the ground state amplitude against the 
number of exponentials N = N in the fits, for the RWcd and point sources. 



(for? =y Q -x Q >0) 

' y,x 



G\ q , (y,x) n\y) 5(A) ( °^ x) ° ) 5(A) Q(x) 



-ik-(y-x) 



(4.9) 



Note that due to the use of mNRQCD for the b quark, the physical momentum p is related to the 
lattice momentum k by 

p = k + Z p ym b v, (4.10) 

where Z p 1 is the renormalisation of the external momentum. Similarly, the physical energy 
Po = Eg of the B meson is related to the energy E v (k) obtained from the fit by 



E B (p)=E v (k)+C v 



(4.11) 



where C v is the velocity-dependent energy shift. Both C v and Z p have been calculated both non- 
perturbatively and perturbatively in Ref. 
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Figure 10: Comparison of relative errors of B meson two point functions from the RWcd, RW and point 
sources, at v = (left) and v = 0.4 (right). 



In order to obtain the random wall correlator, we define 

G£(W)=£ ( G ^' X) ^3(A)***fl(*)^(*)- (4-12) 

We then combine this with the zero-momentum random wall light quark propagator from the same 
^ p to obtain the B meson random- wall correlator: 

1 1 1 " z ~ r ~ i 

C B , RW (a) = -73- £ £ G£ (y,xo,0)* tr ^(y)S(A)G£(y,;to,*) <T^. (4.13) 

y L n z p=\ j l j 

In ( |4.13 ), "tr" denotes a trace over spinor indices only. 

We also compute correlators with gauge-invariant Gaussian smearing for the heavy quark at 
the source and/or sink. The smearing is performed via the operator 



1 + — 4 (2) , (4.14) 
n s 



where is a covariant lattice Laplacian and a, ns are the smearing parameters. The operator 
( 4.14[ ) is inserted in Eq. (4.12) to the left of Q(x) % p (x) (for source smearing) and/or to the left of 
G Vv (y,x) (for sink smearing). 

In Fig. |T(] we compare the relative errors of B meson two-point functions (without smearing) at 
v = and v = 0.4 from the three sources. Similar to vector mesons, an advantage of using random 
wall sources is only seen at time slices less than 10. 

For the heavy-light two-point functions, we use matrix fits with local and smeared sources. 
The fit function has the form 

C/(? ) ^)=£ 1 A^V^ + (-ir 1 £ 1 A^* e -^, (4.15) 



n=0 n=Q 



where the index s (s 1 ) labels the type of smearing at the source (sink). As described in Sec. 4.1 



we 



actually use the logarithms of the energy splittings and the relative excited state amplitudes as the 
fit parameters. 
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0.35 iii 




t Number of non — oscillating states 



Figure 11: Left: A constrained 6+6 exponential matrix fit for the B meson at v = 0.4 (RWcd source; fit 
range? —2. ..32). Right: Fit results for the local B meson amplitude atv = 0.4vs the number of exponentials 
N = N in the Bayesian fit, for the RWcd and point sources. 



Fig. 1 1 shows an example of a 2 x 1 matrix fit of B meson correlators at v = 0.4 with N = N = 6. 
In this fit, comparing the RWcd source with the point source, we see an improvement factor of about 
1 .4 for the ground state energy and local amplitude. Generally the improvement factors we found 
for B mesons are small (no improvement in some cases). 



4.3 Three-point functions 

For xq > yo > z.o, and writing t = xq — yo, T = x$ — zo, the exact all-to-all three-point correlator 



is 



C FJB (t, T, k, p>) = e-^e-^-^e^ Tr 



Gl(y,x)F(x) &(y)f 5 



(4.16) 



where F(x) = 1 for a pseudoscalar meson in the final state, F{x) = (— \) x ' yi for a vector meson 
in the final state, and J? denotes the gamma matrix / derivative operator content of the heavy-light 
current. We now define the sequential-source heavy-quark random wall propagator, based on the 
light spectator quark random wall propagator G p x , (z,xo,—p'): 



G p H (y,zo,x ,k,p')=Z \^ G ^ Z) ° j S(A)f 5 e fk -n(z)G p Xql (z,x () ,-p'). (4.17) 

The random-wall three-point correlator for a pseudoscalar meson in the final state is then 

_^ 111" 2 - ^ r ~ - 

CsjB,JM(t> r > k > P') = -7J- £ I G p q (y,x ,0y -tr Q*(y)% f G^(y,z ,x ,k,p' 
i Li n z „ = i n L 



x e 



p 

-i(k-p')-y 



(4.18) 
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Figure 12: Fits to B — > K three-point functions at 
zero recoil (temporal vector current) from RWcd 
and point sources (data points coincide). 
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Figure 13: Fits to B — > K* three-point functions at 
zero recoil (current: sOojYsb) from RWcd and point 
sources. 



and for a vector meson in the final state we have 

-^111 " z ~ r ~ - -* i 

CjjB,Mh T,k,p>) = -jj-^^G^(y,xojr-tv[rY 5 nHy)f 5 ^G^y,z ,xo,k,p') 

xe -rt-p')-y_ (419) 
For smeared three-point functions, we insert the operator ( |4 .14 ) in the equation for the sequential- 



source heavy-quark propagator ( |4.17D , to the left of Q (z) G^, (z,xq, —p'). 



For the three-point correlators, the fit function has the form 

N F -lN B -\ 
n=0 m=0 



C FJB (p,k,t,T)= £ £ Ahn B) e- Fnt e- Em{ - T -^ + oscillating terms, (4.20) 



and again we actually use the logarithms of the energy splittings and the relative excited state 
amplitudes as the fit parameters. 

Fig. [12| shows fits to B — > K three-point functions at p' = k = 0, v = 0, for the RWcd and point 
sources. These fits have Nb = Nb = Nk = 5 and N% = 0. Correlators with T = 1 1 and T = 12 are 
fitted simultaneously and the range for t is from to T. Comparing the RWcd source to the point 
source, we find an improvement factor of about 1.6 for the amplitude Aqo- 



Fits of B — > K* three-point functions at p' = k = 0, v = 0, are shown in Fig. 13. Here, the 
numbers of exponentials are Nb = Nb =Nk* =Nk* = 6 and the fit range is T = 13, 14; t = 2...{T — 
1). A factor of about 1.3 improvement is observed for the amplitude Aoo for the RWcd source 
compared to the point source (but recall that more inversions were used for the RWcd source). 



5. Simultaneous fits and preliminary form factor results 



The most accurate results for the form factors can be obtained by fitting the two-point and 
three-point functions described in sections 1.1, 4.2 and 03] simultaneously. In the simultaneous 
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Figure 14: Preliminary results for the form fac- 
tors /o, /+ for B — > K decays, obtained from simul- 
taneous non-Bayesian fits with a wide range of T 
in the 3-point function. The left-most points have 
v = (0.4,0,0), k = and p' = 2n/L ■ (-1,0,0). 
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Figure 15: Preliminary results for the form fac- 
tor fj for B —> K decays, obtained from simulta- 
neous non-Bayesian fits with a wide range of T 
in the 3-point function. The left-most points have 
v= (0.4,0,0), k = and f = 2%/L- (-1,0,0). 
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Figure 16: Preliminary results for the form factors T\ , T2 for B — » K* decays, obtained from simultaneous 
Bayesian fits with three values for T. The left-most points have v = 0, k = and p' = 2%/L - (—1,0,0). 



fits, the three-point function ( |4.2Q ) shares the energy parameters F n , F n with the light-meson two- 
point function, and the energy parameters E m , E m with the B-meson two-point function. We fully 
take into account correlations between all data points. 

When multiple values for T in the three-point functions are included, we find that the results 
for the B meson energy and two-point amplitude are significantly more accurate compared to fits of 
the B meson two-point functions alone. This effect is particularly pronounced for the more precise 
B — ^ K three-point functions. 

We have computed the three-point functions for all values of T from to 26 (at v = 0) and to 
20 (at v = 0.4), so that we can investigate which range for T gives the best fits. This investigation 
is still ongoing. Bayesian fits with the full range t = 0...T in the three-point functions turned out to 
be much more difficult once more than 3 or 4 values of T are included (it seems to be impossible 
to achieve # 2 /dof w 1 in this case). 

Therefore, we have also performed non-Bayesian fits with only 1 + 1 exponentials for each 
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meson, including all available values of T and skipping enough points near the sources so that 
the contamination from excited states is seen to be negligible. This only works for B — > K, where 
the signal is still good at large times, and even there only at low recoil. Preliminary results for 
fo, f+ and fx obtained with this method are shown in Figs. |l4| and [l5| respectively. The error 
estimates are from bootstrap. At zero recoil, we get /o(^max) = 0.869(17) from the point source and 
/o(<Zmax) = 0.889(12) from the RWcd source, an improvement by a factor of about 1.4. However, 
at v = (0.4,0,0), k = and p' = 2n/L • (—1,0,0) the point source actually gives more accurate 
results. 

For B — > K* , we used Bayesian fits with three values for T, e.g. T = 13,14,15 and t = 
2 ... (T — 2) for T\. The preliminary results for the form factors T\ and T2 are shown in Fig. 16 We 
have computed the correlation functions for B s — > <p as well, but these still need to be fitted. 

All form factor results presented here include the 1-loop radiative corrections in the heavy- 
light operators as discussed in Sec. ^ a value of a s = 0.3 was used here. The l/m corrections will 
be included once we have the perturbative results for the mixing coefficients. 



6. Discussion 



The stochastic source method we tested reduces statistical errors much more effectively in 
light pseudoscalar mesons than in vector mesons or heavy-light mesons. Its effectiveness is further 
reduced for non-zero momentum, and hence for lower q 2 . For each momentum, and in the case 
of vector mesons also each polarisation, additional inversions are needed for the stochastic source 
method. Our preliminary results with the random wall sources are generally not as as good as in 
[0], which may be due to the lower light-quark mass used here. In calculations with many meson 
momenta, simply increasing the number of point sources may be favoured over using stochastic 
sources if the total computer time is fixed. 
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